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For more than twenty years, the spatial entanglement of photon pairs have been used to uncover interesting 
aspects of quantum physics. A key advantage of spatial variables of single photons relates to the intrinsic high 
dimensionality of its state space, which allows for the investigation of multi-dimensional entanglement and the 
efficient encoding and processing of quantum information. As a drawback, the characterization of the spatial 
correlations over a wide range of spatial modes can be costly, implying a large number of measurements. Here 
we develop and experimentally test novel entanglement criteria based on a periodic discretization of the spatial 
variables. This discretization allows for the definition of d-outcome measurements of arbitrary dimensionality 
that we experimentally implement using a finite set of d periodic spatial masks as mode analysers. Our entan¬ 
glement criteria are computable from 2d^ joint projective measurements, required to characterize the correlation 
of the photons over the spatial masks positioned in the image plane and the far-field of the two photon source. 
Spatial entanglement can then be identified with as few as 8 different measurements. We study the spatial cor¬ 
relations and entanglement detection as a function of the dimensionality d and periodicity of the spatial masks, 
the latter working as a free parameter that can be used to independently tune the measured correlations and 
optimize the entanglement detection. Our results are useful for characterizing spatial entanglement, and have 
applications in quantum cryptography and imaging. 

I. INTRODUCTION 

The transverse position and momentum of single photons 
are utilized as conjugate quantum variables in many interest¬ 
ing fundamental investigations of quantum mechanics such 
as complementarity relations am, continuous variable (CV) 
entanglement Eia and the spatial demonstration of Einstein- 
Rodolsky-Rosen (ERR) correlations ISV ITTIl . In the CV regime, 
the observation of spatial entanglement of photon pairs from 
spontaneous parametric down-conversion (SRDC) typically 
requires the measurement of the detection probabilities asso¬ 
ciated with the joint distribution for the position and momenta 
of the photons, from which non-separability EHU andEPR 
Eiini criteria can be tested. Albeit fundamentally contin¬ 
uous, these distributions in a real-world experiment are sub¬ 
jected to the coarse-graining imposed by the detector reso¬ 
lution lITSl . Even though typical CV entanglement and ERR 
criteria have been adapted for coarse-grained measurements 
|[T4l[T5l, the experimental assessment of the position and mo¬ 
mentum correlations necessary for entanglement detection re¬ 
quire a large number of measurements, as the spatial entan¬ 
glement generated in the SRDC process can be very high¬ 
dimensional. 

Recently, Howland et al csi have shown that the use of 
compressive sensing techniques allows for the reconstruction 
of the SRDC joint detection probabilities with an efficiency 
improvement over a raster scanning procedure with equiv¬ 
alent resolution. This efficiency improvement is twofold: 
firstly, due to the sparsity of the joint detection probabilities 
(the detection modes are highly correlated), a lower number 
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of measurements is needed in the compressive reconstruc¬ 
tion algorithm; secondly, the random binary patterns used in 
the compressed measurement scheme permits the use of a 
higher, approximately uniform joint detection rate. Although 
this scheme allows for a more efficient characterization of the 
spatial correlations, the identification of entanglement is still 
bound to the application of the typical CV entanglement cri¬ 
teria 03 , based on the evaluation of the moments ifTsHSOl or 
entropy EIIES of the reconstructed distributions. 

In the present contribution, we show that a finite set of pe¬ 
riodic amplitude masks (playing the role of mode analyzers) 
can be successfully used to investigate spatial correlations of 
photons in the transverse linear position and momentum do¬ 
main. To this end, we choose a complete set of orthogonal 
mode analyzers comprising d spatial masks with the shape of 
a periodic array of slits, characterized by a slit periodicity T 
and a slit width s = T/d. In this scenario, the spatial vari¬ 
ables of photons are subject to a d-outcome measurement that 
is equivalent to a periodic binning across one of the transverse 
spatial directions. Differently from analysing a d-dimensional 
subset of all available modes ||23l[24]|, the spatial masks act as 
mode analyzers over the entire field of view of the photons. 
Our d-outcome measurement can be regarded as a discrimi¬ 
nation of d-dimensional, orthogonal single-photon images of 
the spatial masks 1^ . 

Based on these d-outcome measurements performed in po¬ 
sition and momentum domains, we derive and implement 
entanglement criteria for bipartite CV systems. We gener¬ 
ate spatially-entangled photon pairs with a SRDC source and 
measure their transverse spatial variables according to this 
discretization procedure. In our experiment, these measure¬ 
ments are realized by positioning the spatial masks in the im¬ 
age plane and far-field of the SRDC source. Our entanglement 
criteria are based on a correlation quantifier - often used as a 
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correlation measure for discrete quantum systems - called the 
mutual predictability |[26l or the coincidence rate IZTl . which 
are directly computable from the 2d^ joint detection rates over 
the spatial masks. Thus, they do not rely on the reconstruction 
of the position and momentum joint probability distributions. 
Interestingly, the periodicity of the spatial masks used in po¬ 
sition (Tx) and momentum (T^) measurements work as free 
parameters that can be independently tuned to optimize the 
entanglement detection. Moreover, we show that useful CV 
entanglement criteria are obtained even for d = 2, which rep¬ 
resents the most extreme coarse-grained measurement possi¬ 
ble with this scheme. It is worth noting that our results are 
applicable to any quantum continuous variable, and can be 
used even if a complete CV measurement has been made, as it 
is possible to post-process the measured distributions into the 
necessary periodic binning structure. 

We tested for spatial entanglement of the photon pairs from 
SPDC using measurement dimensionality ranging from d = 2 
up to d = 20. The data together with the adjusted theoretical 
procedure allow us to reliably confirm CV entanglement with 
as few as 8 joint projective measurements (for d = 2) and 162 
joint projective measurements (for d = 9). In total, we tested 
7344 different combinations of mask periods and dimensions, 
being able to show spatial entanglement in 4432 of them, 60% 
of all the studied cases. 

The paper is structured as follows: in section]^ we provide 
the theoretical background necessary for the development of 
our entanglement criteria; our experimental scheme and mea¬ 
surements with the periodic spatial masks are described in sec¬ 
tion ||n| section ||V] is devoted to the derivation of our CV en¬ 
tanglement criteria; the analysis of our experimental data with 
the derived entanglement criteria is presented in section|V| we 
finally conclude in section [Vl| 


II. THEORETICAL BACKGROUND 
A. The mutual predictability 

The mutual predictability 1^1^ . a building block of our 
entanglement detection procedure, is a correlation function as¬ 
sociated with a measurement of two d-outcome observables 
d and b performed on the two parties of a bipartite quantum 
system Ha C) Hb- Given the joint probability distribution 
Vab{k^ 1) encoding the probability that subsystem A is found 
in the eigenstate \ak) of observable d and subsystem B is 
found in the eigenstate \bi) of observable 6, the mutual pre¬ 
dictability reads 

d-l 

Cab = y2rab{k,g{k)). ( 1 ) 

k=0 

The function ^(-j is a permutation of indices maximizing the 
right hand side of SO that Cab measures the probability 
of obtaining paired (related by / = g{k) (271) measurement- 
outcomes a/c and The mutual predictability attains its 

maximal value equal to 1 for maximally correlated measure¬ 


ment outcomes, whereas its minimal value 1/d is achieved for 
completely uncorrelated outcomes. 

Spengler et al. flSIi showed that entanglement criteria for 
discrete d-dimensional quantum systems can be constructed 
using the mutual predictabilities Q associated with measure¬ 
ments in mutually unbiased bases (MUBs) (23. For this, two 
sets {di} and {bi} of maximally incompatible observables 
(with mutually unbiased eigenbases (291 ) need to be consid¬ 
ered. The entanglement criteria is then built from the linear 
superposition of the mutual predictabilities Ca^bi (^1^ 

m ^ 

= <1 + ^^, ( 2 ) 

i=l 

where “m” is the number of MUBs used. All separable d x d- 
dimensional states satisfy the inequality ([^, and the bound is 
reached by classically correlated states. The use of only two 
MUBs (m = 2) is enough to test for entanglement in all pure 
states, although the use of more MUBs (m > 2) increases 
the distance of Im,d from its upper bound, making more 
robust to noise. In an experimental demonstration using pho¬ 
ton pairs entangled in their orbital angular momentum degree 
of freedom, it was shown that this protocol is able to capture 
bound entanglement (^ . It is worth noting that the entangle¬ 
ment criteria defined in Eq. ^ was generalized in Ref. (^ 
to mutually unbiased measurements (SI. 

In section |IV] we will present a modified version of the in¬ 
equality adapted for the measurements of the spatial vari¬ 
ables of photons with the periodic spatial masks. 

B. Continuous variable entanglement 

It is well known that position \x) and momentum \p) eigen¬ 
states of a CV system satisfy the overlap relation {x\p) = 
e*^^/^/V27rd, and constitute a pair of mutually unbiased 
bases spanning an infinite-dimensional Hilbert space (321. 
Here we consider a bipartite CV state qab with position 
and momentum operators satisfying [xj^pk] = ihdjk (j^k = 
A, B). The joint probability distributions for the position and 
momenta of A and B are given by: 

V{xa,xb) = {xa,xb\Ab\xb,xa), (3a) 

'P{pa,Pb) = {pa,Pb\Ab\pa,Pb)- (3b) 

Given the distributions ( |3a| ) and CV entanglement can be 
identified by applying entanglement criteria based on the mo¬ 
ments (T 8 H 20 I or entropy (2Tl[22l of these distributions. EPR 
steering criteria (T2l[33l[34l are also computable from ( [3a| ) and 
Furthermore, coarse-grained versions of CV entangle¬ 
ment and EPR steering criteria have been recently developed 

Cl El. 

In this work, we consider conjugate spatial variables of sin¬ 
gle photons, namely the transverse position and the wave- 
vector, which satisfy the above commutation relation with 
h = 1. For the sake of brevity, we will refer to the trans¬ 
verse wave-vector as the transverse momentum or simply the 
momentum. Our results, nevertheless, are applicable to any 
quantum continuous variables. 
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FIG. 1: a) Experimental setup. Spatial correlations of entangled photon pairs from SPDC are analyzed with periodic spatial masks displayed 
on SLMs. The blue telescopic imaging system is used for measurements in the position domain, whereas the red Fourier transform lens system 
is used for measurements in the momentum domain. BS: beam splitter; SLM: spatial light modulator; SPAD: single-photon avalanche diode; 
BBO: Beta Barium Borate non-linear crystal, b) Examples of spatial mask sets {Mk} with period T and dimension d = 2 (up) and d = 3 
(down). The spatial masks are defined by a periodic array of slits with slit periodicity T and slit width s — T jd, where the parameter d defines 
the dimensionality of the discretized, d-outcome measurement. 


III. EXPERIMENTAL SCHEME 

In the experimental part of our investigations, we use an 
SPDC source to generate photon pairs entangled in their trans¬ 
verse spatial variables. As mentioned in the introduction, we 
analyze their transverse position and momentum using a finite 
set of orthogonal spatial masks (k = 0,..., d — 1) 

«.(«T) = {;; + (4) 

Each of the mask functions defined in Eq. 0 is a periodic 
square wave with spatial period T and slit width s = T/d. 
The variable z over which the mask function is defined rep¬ 
resents the CV being measured. In our case, it corresponds 
to the horizontal spatial variable at the detection plane: the 
image plane (position) or the far-field (momentum). Note 
also that the particular functions in (|^ describing the distinct 
masks are equivalent upon a suitable translation of coordinates 
in G Z) 

Mk{z-ns;T) = Mk+nimodd){z;T). (5) 

Our experimental setup is sketched in Eig. 0a). We ad¬ 
just the SPDC source to emit frequency-degenerate down- 
converted photons at A = 650nm in a collinear configura¬ 
tion. The down-converted photons are split by a beam split¬ 
ter and directed to spatial light modulators (SLM), which 


are programmed to display the spatial masks described by 
Eq. 0 in the horizontal direction x, as illustrated in Eig. 
0b). A switchable lens system placed before the beam split¬ 
ter is used to produce either the image or the Eourier trans¬ 
form of the SPDC source on to the SLMs. The imaging sys¬ 
tem is characterised by an optical magnification of M = 5, 
while the optical Eourier transform system has an effective 
focal length /e = 333mm. The optical magnifications both 
for image plane (IP) and far-held (EE) configurations were 
chosen in order to maximize the spatial resolution of our 
measurements while keeping the down-converted beams en¬ 
closed within the SLM panel, which have an active area of 
15.36mm x8.64mm. The photons reflected from the SLMs 
are lens-coupled to multi-mode optical flbres connected to 
single-photon avalanche diodes (SPAD). Eurther details about 
our SPDC source and experimental apparatus are provided in 
Sec. I of Supplementary Material. The spatial masks dis¬ 
played on the SLMs work as aperture functions, reflecting the 
photons incident upon the white stripes (where = 1) and 
discarding the photons incident upon the black region (where 
Mk = 0). In Eig. [^b) we illustrate two sets of spatial masks 
with the same periodicity but different dimensionality. 

The effective scaling transformations relating the transverse 
position “x” and momentum “p” coordinates of the source 
with those of the detection plane (denoted further by x^^ and 
x^) are thus of the form x x^^ = Aix and p i-A = ap. 
The parameter a = A/e/(27r) with dimension of length 
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(b) Momentum, 6=2 
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(c) Position, d=4 
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(d) Momentum, d=4 
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FIG. 2: Left: Coincidence counts as a function of the mask displacement ford— 2 and d — A, both for position and momentum measurements. 
In these measurements, the spatial mask #0 is scanned on SLM B across one mask period for each spatial mask displayed on SLM A. Blue 
squares correspond to mask #0, red circles to mask #1, green diamonds to mask #2 and yellow stars to mask #3. From left to right, the 
used mask periods are Tx,Tp = 0.24mm, 0.96mm and 3.12mm. Right: Joint probability distributions over the spatial masks obtained from 
the corresponding coincidence counts on the left. The associated mutual predictabilities are shown below. 


squared is a constant resulting from the action of the optical 
Fourier transform system (see Supplementary Material, Sec. I 
A). The commutation relation for observables associated with 
the rescaled coordinates is: 

= ic^ c = A4a. (6) 

From now on, we will work only with the coordinates of the 
detection plane, and whenever there is no ambiguity we will 
omit the labels “IP” and “FF”. Since the spatial masks op¬ 
erate at the detection planes (never deal with the transverse 
position and momentum per se) we associate the periods T^ 
and Tp with the spatial masks used in IP and FF measure¬ 
ments, respectively. Analogously, we denote Sx = Txld and 


Sp = Tpld the corresponding mask slit widths. Note that both 
Tx and Tp have dimension of length. 

As a first test of the spatial correlations according to our 
measurement scheme, we record the coincidence counts ob¬ 
tained as a function of the relative position of the spatial masks 
on photon A and B. For each mask displayed on SLM A, we 
scan mask #0 on SLM B by the distance of one period, both in 
the image plane and far-field of the SPDC source. These mea¬ 
surements are shown on the left-hand side of Fig. [parts (a) 
and (b) for (i = 2; parts (c) and (d) for d = A]. The strength of 
the spatial correlations is reflected in the coincidence counts, 
whose visibility increases with increasing values of the mask 
periods (from left to right: Tx^Tp = 0.24mm, 0.96mm and 
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3.12mm). For T^^Tp = 0.24mm (equivalent to 30 pixels of 
our SLM), no oscillations were observed. On the other hand, 
for Tx^Tp = 3.12mm, we observe an approximately triangu¬ 
lar oscillation of the coincidence counts, resembling the con¬ 
volution of two slits, each of width Sx or Sp. 

On the right-hand side of Fig. we present the discrete 
joint probability distributions Vxx{k,l) and Vpp{k,l) result¬ 
ing from the measurements with the spatial masks in position 
and momentum domain respectively. These distributions are 
given by the integrals of the joint detection probabilities over 
the spatial masks: 


II 

/ dXA i 

' dXBV\p{xA,XB), 

(7a) 


I ^k-,Tx 

^l]Tx 


Vpp{k,l) = 

/ dxA 1 

' dxB'Pw{xA,XB)- 

(7b) 


I ^k-,Tp d 

^l-,Tp 



The notation M^-t means that the range of integration is re¬ 
stricted by the constraint Mk{x; T) = 1. The joint detection 
probabilities at the SLMs, Vip and Vpp, associated with IP and 
FF measurements, correspond to the rescaled versions of the 
distributions in Eqs. ( [3a| ) and ( [3b| ), respectively. As the pho¬ 
tons at the SPDC source are position-correlated, maximum 
coincidence counts in the image plane occurs for a combi¬ 
nation of masks obeying I = gx{k) = k. In the momen¬ 
tum domain, however, the photons are anti-correlated, so that 
maximum coincidence counts in the far-held are found for a 
combination of masks satisfying I = gp{k) = d — 1 — k. For 
increasing values of the period, the d x d matrices representing 
Vxx{k^ 1) and Vpp{k^ 1) acquire a diagonal and anti-diagonal 
structure, respectively, rehecting the EPR-type spatial corre¬ 
lations displayed by the SPDC photons lO . 

In the next step, we calculate the mutual predictabilities 
Cxx and Cpp corresponding to the discrete joint distributions 
( fTa] ) and ( [75| ), respectively. After substituting in Eq. 0 the 
permutations gx and gp dehned above, we obtain: 


d-1 

'^'Pxx{k,k), 

(8a) 

k=0 


d-1 


'^'Ppp{k,d- 1 - k). 
fe=0 

(8b) 


The experimental values of Cxx and Cpp are given in Fig. 
below the corresponding joint probabilities. For the smallest 
value of the mask period, the measured mutual predictabilities 
are all close to 1/d, indicating no correlations at all between 
the measurement outcomes associated with the different spa¬ 
tial masks. It is also evident from Fig. [^that the correlations 
in the momentum domain are stronger than those in position. 
This happens due to the optical magnihcations used to obtain 
the image and far-held pattern (see Supplementary Material, 
Sec. I A). 


IV. ENTANGLEMENT CRITERIA 

Our goal is to use the mutual predictabilities given in Eqs. 
( [8a| ) and ( [8b| ) to test for entanglement of the spatial variables 


of the two-photon state from SPDC. To this end we need to 
derive the relevant counterpart of the entanglement criteria 
given by 0. The right-hand side of inequality 0 repre¬ 
sents an upper bound for the sum of the mutual predictabil¬ 
ities of d-dimensional, separable states. Derivation of this 
bound is based on the unbiasedness assumption fulhlled by the 
m MUBs at which the joint probabilities are measured 1261 . 
While using two MUBs (m = 2), the bound in inequality Q 
given by 1 + l/d is a function of the dimension parameter 
only. 

Although the position and momentum bases constitute a 
pair of MUBs in the CV scenario 1^ . our d-outcome posi¬ 
tion and momentum measurements do not strictly satisfy the 
condition of maximal incompatibility |[35l for a d-dimensional 
Hilbert space. Then, in order to construct an useful and reli¬ 
able entanglement criterion from the mutual predictabilities 
( [8a| ) and ( [8b| ), we need a tight upper bound for their sum, 
valid for separable CV states. Similar bounds, present in 
other known CV entanglement criteria, are usually obtained 
from uncertainty relations associated with the conjugate vari¬ 
ables of the CV system. For example, in the case of continu¬ 
ous measurements of position and momentum variables, the 
relevant uncertainty relation comes from the Fourier trans¬ 
form relationship between position and momentum eigen¬ 
states (canonical commutation relations). Modifications of the 
Heisenberg uncertainty relation applicable to coarse-grained 
measurements were provided in Ref. (1311361. To deal with 
our d-outcome measurements with the periodic spatial masks, 
we will use the uncertainty relation for the characteristic func¬ 
tions, which has recently been derived in Ref. Ell. 

The operators associated with the measurements involving 
spatial masks can be written as: 

flfc = [ dx^^Mk{x^^\Tx)\x^){x^^\, (9a) 

Jr 

ilk = ldx^Mk{x^-,Tj,)\x^){x^\. (9b) 

dR 


The detection plane bases are related to the posi¬ 
tion/momentum bases as \x^^) = M~^^‘^\x) and \x^^) = 
These two sets of projection operators, {H/c} and 
{U/c}, form complete sets in the meaning that they are suffi¬ 
cient to detect any non-zero position or momentum compo¬ 
nent of the physical state, i.e., Ylk=o 6/^=1 = Ylk=o 

Our CV entanglement criteria, conceptually similar to those 
in Eq. 0, shall be of the form 

Id{Tx^Tp) = Cxx + Cpp ^ Qd{Tx^Tp)^ ( 10 ) 


where the inequality part of (10) is necessarily satisfied when 
qab is separable. A violation of this inequality implies that 
the bipartite CV state is entangled. Our aim is thus to find a 
non-trivial bound 1 + 1/d ^ Qd{Tx,Tp) ^ 2, sharp enough 
to be useful in experimental applications. 

As the starting point we observe that Id is a linear functional 
of the quantum state, so one can restrict further calculation 
of the bound to separable pure states qab = 1^)^ C) |^)^- 
Moreover, the initial steps from the proof of Q, (Eqs. (9) and 
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( 10 ) from ( 261 ), directly apply to the current situation. We 
immediately obtain 

d-l 

C,, + Cpp^J2iPl + 1k)^ ( 11 ) 

k=0 

where pk = (5'| life |4'), qk = (^'1 (^fe 1 ^') and |5') is a single 
photon pure state. 

Theorem 3 from (37) says that 

d—l oo 

E (pI +Ql)^dJ2 \fN\^l3{^7r^N^c/T,Tp), (12) 

k=0 N= — oo 


where 


B(7) = 2V2 


V 2 - yi - cos ( 7 ) 

1 + cos ( 7 ) 


(13) 


and /at are the coefficients of the Fourier series expansion of 
the basic mask function Mq. The constant c comes from the 
commutation relation ([^ between the coordinates of the de¬ 
tection plane. The desired bound Qd{Tx^Tp) is given by the 
right hand side of with the explicit form of the Fourier se¬ 
ries components /at = (1 — where u = 27riN. The 

explicit expression for Qd{Tx^ Tp) is provided in Eq. (SIO) of 
the Supplementary Material. As can be seen, Qd is a function 
of the product of the mask periods Tx and Tp. Interestingly, 
the relation TxTp = ^'Kcjn (n G N) lead to the trivial value 
Qd = 2, meaning that no entanglement detection is possible 
using this specific combination of mask periods. The appre¬ 
ciation of this fact comes from the commutativity of the mea¬ 
surement operators in Eqs. ( |9^ and ( |9b| ), a joint eigenstate of 
which, given by a normalized version of a Dirac comb, is pos¬ 
sible when the periodicity of the spatial masks match those of 
the comb in both domains (TTII^ . 

An infinite sum constituting the upper bound in question 
can be efficiently calculated with the help of a decreasing se¬ 
quence of truncations 


No 

Q^°iT,,Tp) = 2 - d ^ \fMf[2 - B{4Tr^N^c/T,Tp)], 

N=-No 

(14) 

which upper bounds the major quantity Qd{Tx,Tp). The 
derivation of the above truncations is straightforward. Since 
S(-) is never bigger than 2 , one only needs to upper bound 
the terms B{A7v‘^N‘^c/TxTp) in ( [T^ whenever |A^| > A'o- 
The expression ( p^ is recovered with the help of the identity 
Z] at I P = 1 /^’ summation running from — oo to -|-oo. 
Detailed accuracy analysis of the truncations ( p^ is provided 
in Sec. 2 of the Supplementary Material. It is worth notic¬ 
ing that, by using an appropriate number of summation terms 
No, it is always possible to achieve an upper bound which is 
equivalent to Qd {Tx , T^) up to the precision of the experimen¬ 
tal data. 


V. RESULTS 

In plots (a) and (b) of Fig. we show experimental values 
for the mutual predictabilities Cxx ( [ 8 a| ) and Cpp ( [ 8 b| ) obtained 


from measurements using d = 2, 3,4, 5 and 6 . The plot’s 
abscissas are the mask period (left) and slit width (right). In¬ 
terestingly, when plotted against the slit widths, the mutual 
predictabilities display the same behaviour regardless of the 
measurement dimension d. Discrepancies can only be noticed 
for small slit widths, when the mutual predictabilities tend to 
the value of 1/d. This behaviour can be understood by con¬ 
sidering the correlation width of the photons. For a slit width 
s, much larger than the correlation width of the photons, the 
joint probability for the detection of a correlated pair is only 
appreciable within one of the many slits of masks A and B. 
An appreciable joint detection probability in multiple combi¬ 
nation of slits only takes place for mask periods smaller than 
- or on the order of - the correlation width of the photons. For 
instance, considering a Gaussian state and using a slit width 
twice the standard deviation of the joint detection probability, 
approximately 99% of the integrated joint detection probabil¬ 
ity lies within only one combination of slits. Using the param¬ 
eters of our SPDC source and optical setup, we calculate the 
correlation width of the photon pairs in our experiment to be 
^ip{xa — xb) = 0.068mm for measurements in the image 
plane and Aff{xa T xb) = 0.037mm for measurements in 
far-field (see Supplementary Material, Sec. IB). 

For each measurement dimension, we can use any combi¬ 
nation of mask periods in position and momentum domain to 
evaluate Id{Tx^Tp) and test for CV entanglement according 
to our entanglement criteria given in inequality ( p^ . In the 
colormaps (c) and (d) of Fig. we display the experimental 
values obtained for Id{Tx,Tp) as well as the corresponding 
values of Qd{Tx^Tp) for d = 2 and d = 6, respectively. For 
the calculation of Qd, we summed up A'o = 6000 terms in 
the series ( p^ . In that way, we obtained a precision (in de¬ 
termination of the upper bound relevant for the criteria) up to 
the fourth decimal place, which is two orders of magnitude 
smaller than our experimental errors. Our experimental errors 
for Id up to d = Q are of the order of 10 “^, and are calcu¬ 
lated assuming Poissonian statistics for the measured coinci¬ 
dence counts. The dependence of the bound on the product of 
the spatial periods Tx and Tp is clearly seen in the symmet¬ 
ric hyperbolic structures of the colormaps representing Qd. 
Higher-resolution colormaps of these bounds, where the os¬ 
cillating structures are better visualized, are provided in Fig. 
SI of the Supplementary Material. In the right-hand side of 
Figs, [^c) andj^d), we show colormap plots of the difference 
'Id — Qd\ where positive values (in blue) indicate entangle¬ 
ment detection. For small values of the mask periods, no en¬ 
tanglement detection is achieved, as in this case very weak 
correlations are observed both in position and momentum do¬ 
main [see plots (a) and (b) of Fig. [^. For larger periods, 
the observed correlations over the spatial masks are stronger, 
leading to higher values for the mutual predictabilities that al¬ 
low for entanglement detection. Note that, as the correlation 
width of the photon pairs is smaller in the far-field than in the 
image plane configuration, Id{Tx,Tp) grows faster with Tp 
than with Tx. 

In order to study the influence of the measurement dimen¬ 
sion on the detection of entanglement according to inequal¬ 
ity ( p^ , we use the mutual predictabilities given in Fig. 
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FIG. 3: Plots of the experimental values obtained for the mutual predictabilities in (a) position and (b) momentum domain using measurement 
dimensions d = 2, 3,4, 5 and 6. The mutual predictabilities are plotted against the mask periods Tx or Tp (left) and slit width Sx = Txjd or 
Sp — Tpjd (right). In colormaps (c) and (d), we display the experimental values obtained for Id {Tx , Tp) (left), the bound functions Qd {Tx , Tp) 
(middle), and the difference between them (right), where positive values leading to entanglement detection are color mapped in blue. Data in 
(c) corresponds to d = 2 and data in (d) to d = 6. 


to calculate Id{T,T), with T^ = Tp = T, as shown in the 
colormap (a) of Fig. In this graph, the columns corre¬ 
sponding to d = 2 and (i = 6 are equivalent to the diagonals 
of the right-hand side colormaps (c) and (d) of Fig. As 
expected, the smallest mask period for which entanglement 
is detected grows with the measurement dimension. This be¬ 
haviour can be further examined in colormap (b) of Fig. 
where we show a new data set for measurement dimension up 
to (i = 20. Due to the statistical errors of our measurements, 
the largest measurement dimension for which entanglement 
detection was successful is d = 9 (see Table [I]). 

In our experiment, we work with mask periods ranging 
from Tmin = 0.24mm up to T^ax = 8.64mm (being the 


largest mask period equivalent to the vertical length of our 
SLM panel). The strength of the measured correlations, how¬ 
ever, is intimately related to the slit width of the spatial masks. 
This relation is better visualized in plot (c) of Fig. where 
the measurements shown in colormaps (a) and (b) of Fig. 
are plotted against the common slit width Sx = Sp = T/d. 
As before, points for which 'Id — Qd' is positive are de¬ 
picted in blue. In the plot discussed here, an approximately 
measurement-dimensionality independent trend is observed, 
suggesting a cutoff for the mask slit width below which entan¬ 
glement detection is not possible. The two vertical red lines 
indicate the maximum slit width, s^ax = ^max/<^, used in the 
measurements with d = 12 and d = 29. For d = 29, Smax 
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FIG. 4: (a) and (b) Colormaps of the experimental values obtained for ''Id — Qd' as a function of the mask period and measurement dimension, 
using Tx — Tp — T. Spatial entanglement is detected in the blue region where Id > Qd- Colormap (a) corresponds to the data presented 
in Fig. 1^ (c) Plot of the measurements corresponding to colormaps (a) and (b) as a function of the associated slit width s — T/d. (d) 
Each black dot corresponds to a distinct experimental evaluation of our entanglement criteria |T0] ), plotted against the dimensionless parameter 
TxTpj (d^ J\Aa). In total, there are 7344 measurements, 4432 of which lead to entanglement detection. The red horizontal line corresponds to 
7cz — Sd = 0, above which entanglement is detected. 


is clearly smaller than the minimal slit width required to mea¬ 
sure correlations strong enough for entanglement detection. 
Larger slit widths ford = 20 would be achievable for mask 
periods greater than the beam diameter, implying a practically 
null detection probability over some of the spatial masks. We 
notice that in the limit of mask periods much larger than the 
beam diameter, a (i-outcome spatial mask analyzer is in prac¬ 
tice reduced to a 2-outcome measurement equivalent to the 
sign operator — with mask #0 and ^{d — 1) covering the 
positive and negative part of the x axis respectively. 

By analysing our full experimental data in every possi¬ 
ble combination of mask periods and dimension d, we ob¬ 
tained 7344 different measured values of Id{Tx,Tp). In 
plot (d) of Fig. 1^ we display all these values, subtracted 
by the corresponding bounds Qd{Tx,Tp), as a function of 
TxTpI{d‘^Ma). Beside statistical errors for some of the 
points lying close to the red horizontal line indicating Id — 
Qd = 0, we detect spatial entanglement in 4432 out of the 
7344 measurements displayed in this plot. With the purpose 


of providing numerical examples of our experimental results, 
we show in Tablethe largest values of Id{Tx^Tp) obtained 
for each measurement dimension. As mentioned before, vi¬ 
olations of inequality are obtained for measurement di¬ 
mensions up to (i = 9 and some violations smaller than the 
statistical errors are obtained ford = 10 and d = 12. No 
entanglement detection is achieved for d=lbor d = 20. 

It is worth mentioning that s'^ = Txl{Md) and s'^ = 
Tp/{ad) are the slit widths of the spatial mask analyzers 
associated with the position and momentum distributions of 
the photons at the SPDC source. A maximal distance of 
Id{Tx,Tp) from the bound is achieved for a certain value of 
s'^Sp, which we estimate from plot (d) of Fig. 01 to be ap¬ 
proximately s'^Sp ^ 5. Thus, the measurement mstribution 
in plot (d) of Fig. [^indicates the best range of spatial mask 
parameters useful in the application of the entanglement cri¬ 
teria ([To|.. Considering spatial masks in the position and mo¬ 
mentum domain with the same slit width, this maximal dis¬ 
tance corresponds to an optimal slit width of 5opt ~ 0.93mm, 
















































9 


d Ta;(mm) Tp(mm) Id{Ta,,Tp) Qd{Ta,,Tp) 


2 

3.36 

1.20 

1.73 

zb 

0.01 

1.570 

3 

3.36 

2.16 

1.69 

zb 

0.02 

1.525 

4 

7.44 

2.16 

1.72 

zb 

0.02 

1.557 

5 

6.00 

2.64 

1.66 

zb 

0.02 

1.488 

6 

7.92 

4.08 

1.68 

zb 

0.02 

1.546 

8 

7.92 

5.76 

1.61 

zb 

0.04 

1.517 

9 

8.64 

5.76 

1.56 

zb 

0.05 

1.491 

10 

7.20 

5.04 

1.47 

zb 

0.05 

1.438 

12 

7.20 

8.64 

1.49 

zb 

0.06 

1.454 

15 

8.64 

7.92 

1.31 

zb 

0.11 

1.434 

20 

8.64 

8.64 

1.16 

zb 

0.14 

1.386 


TABLE I: Table with the highest values of Id{Tx, Tp) obtained for 
each measurement dimension. The experimental errors are calcu¬ 
lated assuming Poissonian statistics for the measured coincidence 
counts. 


which visually agrees with the peak in plot (c) of Fig. As 
the beam diameter both in position and momentum measure¬ 
ments was chosen to be encompassed within the SLM panel, 
t^max = ^max/-^opt ~ 9-3 gives the largest measurement dimen¬ 
sion that can be used to test for entanglement in our system 
with the optimum slit width Sopt- 

VI. CONCLUSION 

Spatial degrees of freedom of photons are an attractive 
choice for the experimental study of entanglement in high¬ 
dimensional systems. In this contribution, we showed that 
spatial entanglement can be detected from the correlations 


of the photons over a finite set of orthogonal, periodic im¬ 
ages, exploring the full-field of view of the transverse vari¬ 
ables of the photons. Thus, despite the infinite dimensional¬ 
ity of the CV quantum system, our entanglement criteria is 
evaluated from the outcomes of 2d‘^ joint projective measure¬ 
ments, which can be made as few as 8 by using d = 2. These 
image-based entanglement criteria can be understood as a dis¬ 
crimination of (i-dimensional, orthogonal single-photon im¬ 
ages, performed in conjugated transverse variables of the pho¬ 
tons, and may have applications in ghost imaging 13^ and 
CV quantum key distribution 1401 . Our entanglement crite¬ 
ria are applicable to any quantum continuous variable, and we 
believe that they shall be useful for entanglement detection in 
other experimental implementations of CV entanglement EJ- 

m. 
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Supplementary Material for ^^Testing for entanglement with periodic images’’ 

I. SPDC SOURCE AND EXPERIMENTAL APPARATUS 

We use a type-I SPDC source based on a 2mm long BBO crystal and adjusted to emit frequency-degenerate, down-converted 
photons in a collinear configuration. Our pump is a ^ 30mW, vertically polarised He-Cd laser (Kimmon Koha) operating in the 
continuous wave regime at the wavelength of Xp = 325nm. A dichroic mirror is used after the BBO crystal to refiect the pump 
laser whilst transmitting the lower frequency, down-converted photons. Spatially entangled down-converted photon pairs at the 
central wavelength of A = 2Xp = 650nm are selected by the use of lOnm wide interference filters positioned in front of the 
detection system. 

A 50:50 beam splitter (BS) positioned at 400mm from the SPDC source is used to probabilistically split the down-converted 
photon pairs into propagation paths labeled A (reflected output port) and B (transmitted output port). In each propagation path, 
we use phase only spatial light modulators (SLM) to generate the spatial masks used in the measurement of photon A and B. Our 
SLMs (Holoeye PLUTO) panel have an active area of 15.36mm x 8.34mm containing an array of 1920 x 1080 pixels of area 
8/im X 8 /im. The phase pattern displayed on the SLMs is a hologram of the spatial masks given in Eq. (4) of the manuscript. 
The hologram’s diffraction grating is blazed to maximise the diffraction intensity into its first diffraction order, which is coupled 
into multi-mode optical fibres connected to single-photon avalanche diodes (SPAD). This detection system is then equivalent to 
a mode analyser comprised of an aperture function with the shape of the spatial masks Mk . 


A. Imaging and Optical Fourier transform system 

We use an interchangeable lens system to produce either the image or the Fourier transform pattern of the SPDC source onto 
the SLMs. This interchangeable lens system is positioned between the SPDC source and the BS used to split the down-converted 
photons, and thus acts both on paths A and B. Mounted on flip mounts, the optical Fourier transform lens system consists of 
a 200mm focal length lens located halfway from the SPDC source to the BS and the telescopic imaging system comprises a 
50mm and a 150mm focal length lenses. These interchangeable systems are used to Fourier transform or image the generated 
two-photon field onto the BS, respectively. On paths A and B, the SLMs are positioned at an imaging plane of the BS (800mm 
from the BS) by means of a telescopic imaging system comprising lenses of focal length 150mm and 250mm. Overall, from the 
SPDC source to the planes of the SLMs, we have an imaging system with magnification Af = 5 or an optical Fourier transform 
lens system with effective focal length /e = 333mm. 


B. Two-photon state 

With the aim of estimating the transverse correlation width of the photon pairs in our experiment, we give a brief description 
of its spatial structure. The two-photon field from SPDC presents a high degree of correlation between its spatial and temporal 
modes mm. When working in the quasi-monochromatic approximation, a pair of correlated frequency modes is selected by 
narrow spectral filtering of the two-photon field. In this case, the spatial structure of the two-photon field is described by a 
pure state entangled over a large number of spatial modes. Considering a thin non-linear crystal and working in the paraxial 
approximation 13, the transverse spatial state of the generated two-photon field in position representation can be written as m 

\'^)AB=j J '^{PA,PB)\pA)\pB)dpAdpB, (SI) 

where is the transverse coordinate of photon i. The two-photon transverse probability amplitude 4^(py^, p^) is non- 

separable in the transverse coordinates of the photons, and it is approximately factorable in the x and y transverse coordinates: 

'^{pa^Pb) ~ i’{xA,XB)'ip{yA,yB). (S2) 

We are interested in the spatial correlations as probed by the spatial mask analysers displayed on the SLMs. As the spatial masks 
are defined in the horizontal direction x, we consider only the horizontal part of two-photon probability amplitude in Eq. ( [S^ , 
which can be written as n 


'tp{xA,XB) =W{xA+XB)r{xA- Xb), (S3) 

where W{xp) cx is the transverse field distribution of our pump beam (in the x direction) at the crystal, Wp = 0.464mm 

is the pump beam radius and F is the Fourier transform of the phase-matching function which, for collinear phase-matching. 


2 


assumes the form nn 


r(pA,PB) OC sine 


L\p 

Sir 


iPA-Psf , 


(S4) 


where L is the length of the non-linear crystal, and pi is the transverse momentum variable of the down-converted photon i. 
Notice that we are referring to the transverse wave-vector as transverse momentum, which is equivalent to using h= 1. 

The correlation width of the two-photon field on the plane of the SLMs depend on the optical magnifications associated with 
our imaging and Fourier transform system. Let (image plane) and (far-field) be the horizontal coordinates 

at the detection planes of the SLMs A and B, respectively. Considering our imaging system with magnification M = 5, the 
propagated two-photon transverse probability amplitude in the imaging plane is 

(S5) 


For our optical Fourier transform system with effective focal length /e = 333mm, the propagated two-photon probability 
amplitude in the far-field is 


(S 6 ) 

^ ^ \ a a J 

where a = A/e/( 27 r), A is the wavelength of the down-converted photons and is the Fourier transform of t/;, 

^ (S7) 

yj 2 'Ka J 

In our experimental setup, ol ^ 0.0345mm^ maps the momentum variable into position space, = ap, and is associated with 
the optical Fourier transform system. 

As usual, we define the correlation width of the down-converted photons as the standard deviations of their relative position 
(imaging plane) or total momenta (far-field) O. These standard deviations are calculated from the two-photon probability 
amplitude in the imaging plane, Eq. ( [S5] ), and far-field, Eq. ( [S 6 | ). For the calculation of the correlation width in the imaging 
plane, we use the standard deviation based on the sine phase-matching function, Eq. ( [S4| ), obtaining Q 

^ipixA - xb) = 0.068min, (S 8 ) 


whereas the correlation width in the far-field is 

1 OL 

Aff(xa + Xb) =-~ 0.037mm. 

2 Wp 


(S9) 


II. ANALYSIS OF THE BOUND Qd (T,, Tp) 

The upper bound present in the entanglement criteria, given in a more explicit way reads 

r 'i - ^ ^ 1 - cos{2'KN/d) \/2 - yi - cos(27rAr2/7) 

■ddK X. p) - ^2 1 + cos(27rA^V7) 


(SIO) 


where 7 is given by 


T^Tp 

7 = 


27rA^Qf 

The derived truncations forming a decreasing sequence of upper bounds is then given by the expression: 


o (T O^0(T \ ^ ^ ~ cos{2TrN/d) f V2- - cos{2TrN^ /-f) 1 ' 

^d[lx,-Lp)^^d ^2 y 1 + cos(27riV2/7) 


(Sll) 


(S12) 


In the data analysis described in the main part of the paper, all values of Qd were obtained from Eq. ( |S12| ) using Nt^ = 5000. 
In order to provide examples of the fine structure of the function Qd and its dependence on the periods and dimension of the 
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FIG. SI: Colormaps (a) to (d): High-resolution 2.83 Mega-pixels) bounds for entanglement detection calculated according to Eq. ( |S12| ) 
using iVo = 5000 for (i = 2 (a and b) and d = 6 (c and d). (a) Q 2 {Tx,Tp)-, (b) Q 2 {Tx, aTp^)-, (c) Q6{Tx,Tp); (d) Q6{Tx,aTp^). (e) 
Plot of Q 2 (0.96mm, 3.12mm) calculated from Eq. ( |S12| ) (upper dotted line) and Eq. ( |S10| ) (lower dotted line) as a function of the number of 
summation terms Nq. The shaded area correspond to the possible values of the bound, (f) Comparison of the convergence of the bound with 
A^o for mask dimensions c/ = 2, 6 and 20. 


spatial masks, we calculate high-resolution versions of the functions Q 2 {Tx^Tp) and Q^iT^^Tp) displayed in Fig. (4) of the 
main manuscript. These high-resolution versions are shown in Fig. S[^(a), for d = 2, and Fig. 8|^(c), for d = Q. In these 
colormaps, we show Qd with a resolution of 0.005mm, which is 48 x higher than the resolution of our measurements. The 
dependence of Qd on the product of the mask periods [the parameter Eq. ( |S11| )] becomes evident in Fig. 8(^1 (b) and S[^(d), 


where we show, respectively, Q 2 and Qo as a function of and aT^ In these figures, the typical hyperbolic structures that 
appear in Qd become linear patterns. All colormap plots of Id and Qd in this work were made according to Cube Helix colormap 
|6l , having a perception of brightness that increases monotonically with the signal intensity. 

The convergence of Qd with the number of summation terms in Eqs. ( |S10| ) and ( |S12| ) is illustrated in Figs. 8|M](e) and 8 |^ 
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(f). In Fig. S|^(e), we show a plot of Q 2 for = 0.96mm and Tp = 3.12mm as a function of the number os summation terms 
Nq. The lower and upper dotted lines correspond to Q 2 calculated from Eqs. ( |S1Q| ) and ( |S12| ), respectively, and the shaded area 
between them represent the possible values for Q2(0-96mm, 3.12mm) as a function of Nq. After summing up 5000 terms, the 
uncertainty in its determination is of the order of 10“"^. As our experimental errors for Id are of the order of 10“^, we round the 
obtained value for the bound to the third decimal place, giving Q 2 (0.96mm, 3.12mm) = 1.546. In Fig. SSI (f) we compare the 
convergence in the calculation of Qd for mask dimensions d = 2, 6 and 20. Using A'o = 5000, we obtain an uncertainty in the 
determination of Qd smaller than (at least) 10“^ for all used mask dimensions. 
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